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In this work we generalize and subsequently apply the Effective Field Renormalization 
Group technique to the problem of ferro- and antiferromagnetically coupled Ising spins 
with local anisotropy axes in geometrically frustrated geometries (kagome and pyrochlore 
lattices). In this framework, we calculate the various ground states of these systems and 
the corresponding critical points. Excellent agreement is found with exact and Monte Carlo 
results. The effects of frustration are discussed. As pointed out by other authors, it turns out 
that the spin ice model can be exactly mapped to the standard Ising model but with effective 
J3 ' interactions of the opposite sign to those in the original Hamiltonian. Therefore, the ferro- 

magnetic spin ice is frustrated, and does not order. Antiferromagnetic spin ice (in both 2 and 
. 3 dimensions), is found to undergo a transition to a long range ordered state. The thermal 

^ I and magnetic critical exponents for this transition are calculated. It is found that the thermal 

Q I exponent is that of the Ising universality class, whereas the magnetic critical exponent is 

different, as expected from the fact that the Zeeman term has a different symmetry in these 
CN . systems. In addition, the recently introduced Generalized Constant Coupling method is also 

1^ . applied to the calculation of the critical points and ground state configurations. Again, a 

I very good agreement is found with both exact, Monte Carlo, and renormalization group cal- 

culations for the critical points. Incidentally, we show that the generalized constant coupling 
O ' approach can be regarded as the lowest order limit of the EFRG technique, in which corre- 

lations outside a frustrated unit are neglected, and scaling is substituted by strict equality of 
the thermodynamic quantities. 

'a _ 

^ ■ PACS numbers: 75.10.Hk, 75.25.+Z, 75.40.Cx 

O : I. INTRODUCTION 

O 

Geometrically frustrated antiferromagnets (GFAF) have emerged in the last years as a new class 
^ . of magnetic materials with many unusual physical propertiesliliiH. In these materials, the elemen- 
tary magnetic unit is the triangle, which makes it impossible to satisfy all the antiferromagnetic 
bonds at the same time, with the result of a macroscopically degenerate ground state. Examples 
of GFAF are the pyrochlore and the kagome lattices. In the former, the magnetic ions occupy the 
comers of a 3D arrangement of comer sharing tetrahedra; in the later, the magnetic ions occupy 
the comers of a 2D arrangement of comer sharing triangles (see Fig. [1]). In the case of materials 
which crystallize in the pyrochlore structure, the static magnetic susceptibility follows the Curie- 
Weiss law down to temperatures well below the Curie-Weiss temperature. At this point, usually 
one to two orders of magnitude smaller than the Neel point predicted by the standard mean field 
(MF) theory, some systems exhibit some kind of long range order (LRO), whereas others show a 
transition to a spin glass state (SG). This is a striking feature for a system with only a marginal 
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amount of disorder. Finally, there are some pyrochlores which do not exhibit any form of order 
whatsoever, and are usually regarded as spin liquids. In the case of the kagome lattice, even though 
there are very few real systems where this structure is realized, the magnetic properties fall in two 
major categories: the vast majority of the compounds studied show a transition to a LRO state with 
a non coUinear configuration of spins, ailda few systems exhibit no LRO, but a SG like transition. 

Surprisingly, Harris and coworkerstlBLl recently showed that geometrical frustration can also 
arise from ferromagnetic interactions. This cannot occur for Heisenberg spins but, in some cases, 
a strong single ion anisotropy along high symmetry directions of the crystallographic structure can 
force the spins to point in such directions, so they can effectively be considered as Ising spinsB. 
In this case, the Ising spins with ferromagnetic interactions do not order at any finite temperature. 
They analyzed this problem in the context of Ising spins in the comers of a pyrochlore lattice, 
and termed this system spin ice pyrochlore, because it can be mapped to the well know proton 
ordering problem in the ice lattice. Furthermore, in full equivalence with the ice problem, the spin 
ice model presents a slow dynamic, due to the existence of large energy barriers associated with 
rearrangements of the spins in the ground state, resulting in a spin freezing which resembles a spin 
glass, but occurring in a system where chemical disorder is absent. 

While at first sight this can look like a very interesting academic problem, it turned out to be 
also of enormous practical interest, because various pyrochlore oxides were found to have strong 
single ion anisotropy along the (1 1 l)-directions of the lattice. This effect usually occurs for rare 
earth magnetic ions with a large orbital moment, as it is the case of the pyrochlores Ho2Ti207, 
Dy2Ti207, and Yb2Ti207, among others, which are now known to be very good realizations of the 
spin ice model. 

From the theoretical point of view, apart from the interest of the spin ice system itself, this 
system also offers a unique opportunity for testing models which incorporate geometrical frustra- 
tion, for the obvious reason that the Ising Hamiltonian is simpler to study mathematically than 
the Heisenberg one and, in some cases, analytical results are easily obtained, which can be read- 
ily compared with experimental data or simulations. If such a comparison is adequate, one can 
get some confidence that these techniques can also be applied to the Heisenberg problem in the 
pyrochlore lattice, which is an even more interesting problem. 

This is precisely the purpose of the present work. We will study the critical properties of 
geometrically frustrated Is'mgferro- and antiferromagnetshy introducing a generalization of a well 
known real space renormalization group technique, the so called Effective Field Renormalization 
Group (EFRG) method§0, which has been successfully applied to a large variety of systems (see 
Ref. |12| and references therein), and it is known to provide very good estimates of the critical 
temperature and the critical exponents. We apply this technique to two model systems: the spin 
ice problem in the pyrochlore lattice, and a 2D analog to the spin ice problem, which consists of 
Ising spins in the kagome lattice with local axes of anisotropy along the heights of the triangle. 
This second example, though academic, is a very instructive one, as it will serve as an example to 
develop the present method. Our first task will be to generalize the EFRG technique to include the 
effect of geometrical frustration, that is, the non-Bravais character of the lattices considered. We 
will then study the existence of critical points and their values for both ferro- and antiferromagnetic 
interactions. The obtained values are compared with exact values known for the Ising model in the 
kagome latticeliJ and with numerical estimations of this quantity obtained by Harris and Bramwell 
for the spin ice pyrochloreQ. We find that, even at the lowest order of approximation, in which 
correlations with ions further that NN are neglected, the estimates are excellent for both lattices. 
Furthermore, we show that these spin ice-like Ising models can be related to the standard Ising 
modelBil. Therefore, we are able to understand why the ground state given by the ice rules is never 
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(a) Kagome lattice. 



(b) Pyrochlore lattice. 



FIG. 1 : The structures considered in this work. The black thick lines indicate the axes of anisotropy. Spin 
labels are of the form ia, where a is the label of the sublattice and i the index of the spin inside the 
corresponding sublattice. The normal vectors na have been also depicted. 

found in these systems, in agreement with numerical evidences provided by other authorsii. Next, 
we evaluate the thermal and magnetic exponents for the transition to an ordered state that is found 
for antiferromagnetic spin ice, namely, a ground state in which all the spins point out or all the spins 
point in. These exponents are calculated for two limiting cases: the one corresponding to Ising 
spins with a single anisotropy axis, that is, the standard Ising model, and the one corresponding to 
the spin problem. For the standard Ising model the calculated exponents are found to be in good 
agreement when compared with exact values for 2D lattices, and with Monte Carlo simulations 
for the 3D lattice, for the Ising universality class. For the spin ice case, the thermal exponent is 
the same as in the standard Ising model, whereas the magnetic one is found to be different for 
both lattices, in agreement with the fact that the Zeeman term in this case does not have the same 
symmetry as in the standard Ising model. 

Additionally, we find interesting to compare the results obtained from this new method with 
the ones predicted bvthe Generalized Constant Coupling (GCC) model, recently developed by 
the present authorsoEl, which has been shown to provide very accurate results for the magnetic 
properties of Heisenberg GFAF. The results for the critical points and ground state configurations 
obtained in the GCC framework are found to be in excellent agreement with the ones found by us- 
ing the EFRG technique, in spite of the mathematical simplicity of the GCC method. Incidentally, 
we show that this agreement is not a coincidence, but it is related to the fact that the GCC method 
can be considered as the zero-th order approximation of the EFRG scheme. 

The remainder of this paper is organized as follows: in the next section, we introduce the sys- 
tems we will deal with and a brief summary of the EFRG method, together with the corresponding 
generalization to deal with geometrically frustrated lattices. In the third section we develop the 
EFRG scheme for both the kagome and pyrochlore lattice with anisotropic Ising interactions. In 



section |V|, the finite scaling hypothesis is used to calculate the non-trivial fixed points of the RG 
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transformation (the critical points), together with the possible ground state configurations of the 
system. In section ^ we calculate the thermal and magnetic critical exponents for the long range 
ordered state found in the previous section for antiferromagnetic interactions. In section 0, we 
apply the GCC model to the present problem and, again, the critical points and ground state con- 
figurations are calculated and compared with the EFRG method. In section |VII| , the connection 
between both schemes is furnished. Finally, section is devoted to the conclusions. 



II. STATEMENT OF THE PROBLEM 

In the systems we will consider below, the strong single-ion anisotropy forces the spins in the 
lattice to point in high symmetry directions and, to a first approximation, they can be regarded as 
Ising spinsB. For example, in the spin ice pyrochlore, each ion in the elementary tetrahedral unit 
points in a different (1 1 l)-type direction, which connects the spin with the center of the tetrahedra. 
Therefore, the spin variables are not scalar ones, but have vectorial character. We will represent 
each of the spins by a vector = SiJia, where Sia takes the values ±1 and ria stands for the 
unitary vector in one of the directions of anisotropy, where a takes the values a = A, B,C, . . ., 
depending on the number of different directions of anisotropy (see Fig. [I[). The index i runs over 
all the spins with a given normal vector. 

The Ising Hamiltonian with only NN interactions in adimensional units can be then put as 

n = K^^Sia- Sjf3, (1) 

where K = PJ is the adimensional coupling constant, which is positive for ferromagnetic interac- 
tions and negative for antiferromagnetic interactions. We can see that, in this way, we have defined 
a finite number of sublattices, characterized by the corresponding orientation of the Ising spins, 
which has the advantage that a spin in sublattice a interacts only with spins in different sublattices. 

The geometries we will consider in this work are those depicted in Fig. [l]: anisotropic Ising 
spins in the kagome lattice, which can be considered as a 2D analog to the spin ice problem, and 
the spin ice problem, that is, anisotropic Ising spins in a pyrochlore lattice. 

In order to study the critical properties of these systems, we will make use of the so called 
EFRG technique, which has been proved to be a very good RG scheme for dealing with critical 
phenomena. We will not present the method here, as there is an excellent review on this subjectEl, 
but only the most important notions. 

According to the idea of scaling)ll, close to the transition the singular part of the free energy 
scales as 

fsie,H) = l-''Ml'/''ejy"H), (2) 

where e = K — Kc, with Kc the critical coupling, / is an arbitrary scaling factor, d is the di- 
mensionality of the lattice, yn is the magnetic critical exponent and H the applied magnetic field, 
and u is the correlation length critical exponent. Any other thermodynamic quantity, P, can be 
evaluated as derivatives of this singular part of the free energy near the transition, and will behave 
as a power law 
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where cr is the critical exponent of the quantity P. For example, the magnetization, near criticality, 
behaves as 

m(e, H) = l-'^+y"m{l^/''e, F"H). (4) 

On the other hand, according to the finite size scaling hypothesis, the generalized scaling rela- 
tion obeyed bx any thermodynamic quantity P, taking into account the finite size L of the system, 
can be put asll3 

P(e, H, L) = l^Pil^'^'t, l^^H, r^L), (5) 

where (f) = a/u is the anomalous dimension of the quantity P. By choosing a smaller system of 
size L', fixing the scaling factor to the value / = jj, and introducing the variables e' = l^^^e = 
K' — Kc and H' = H, we can rewrite relation (Q) in the form 

Pu{K\H') ^ Pl{K,H) 

This is the basic equation from where a mapping {K, H) {K\ H') can be established, and the 
corresponding critical points and exponents obtained. 

By applying relation @ to different thermodynamic quantities, different RG approaches are 
generated. For example, on applying (||) to the correlation length, we have 

i'{K',H') ^ aK,H) 
L> L ' 

which means that, for the correlation length, the anomalous dimension is the unity. 

The EFRG approach corresponds to apply relation to the order parameter (magnetization). 
The orderparameter for two clusters of different size is calculated by using the Callen-Suzuki 
identityBtl 



where the partial trace is taken over the set of p spin variables specified by the finite size cluster 
Hamiltonian TY^; Op is the corresponding order parameter, and ( . . . )^ indicates the usual canonical 
thermal average over the ensemble defined by the complete Hamiltonian 7i. The idea is then to 
replace the effect of the total Hamiltonian not included in the partial trace in (|]) by fixed values 
of the magnetization of the ions outside the cluster, which act as a symmetry breaking field, h. In 
this way, the order parameter in the finite cluster can be computed for 6 -C 1 and H <^1, and an 
equation of the form 

mp(ir, 6, H) = (Op) = U{K) b + gp{K) H, (9) 

is obtained. Repeating the same calculation for a cluster of different size, p', and assuming that 
both magnetizations are related by the scaling relation (Q) near criticality, we obtain the following 
equation 



fpiK') b' + gp,{K') H' = I'-y" fp{K) b + f-y» gp{K) H. 



(10) 




(a) 1-spin cluster. 



(b) 3-spins cluster. 



FIG. 2: Clusters used in this work for the kagome lattice. Dark spheres represent spins belonging to the 
cluster itself and white spheres the corresponding surrounding sites (those creating an effective field). 



As h and h' are, in some sense, magnetizations, they are also assumed to satisfy @, so we have the 
following relation 



U,{K') = U{K), 



(11) 



from which the fixed points of this RG transformation and, thus the critical points, are found by 
solving this equation for K' = K = K^,. The corresponding thermal eigenvalue is obtained as the 
corresponding eigenvalue of this recursion relation, linearized near the critical point 
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and by making use of ([7|) together with the definition At = / 
correlation length critical exponent as 



V 



1 _ In/ 

The magnetic critical exponent is calculated from 

together with the scaling relation for the magnetic fields, W = H, as 



(12) 

l^/", we can calculate the 
(13) 



Vh 



d 



1 ^^ 9p{Kc, 
In/ gp'{,Kc 



(14) 



(15) 



The only source of inaccuracy in these relations is in the finite sizes of the clusters. However, it has 
been found that these kinds of approaches are applicable to rather small systems and. sometimes, 
they are even more accurate than other RG approaches, for clusters of adequate sizeli3. 

Even though the EFRG scheme has been successfully applied to a broad range of problems, 
the original formulation of this method is not well suited to deal with the effects of geometrical 
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(a) 1-spin cluster. 



(b) 4-spins cluster. 



FIG. 3: Clusters used in this work for the pyrochlore lattice. Dark spheres represent spins belonging to the 
cluster itself and white spheres the corresponding surrounding sites (those creating an effective field). 

frustration, that is, magnetic systems with a non-Bravais lattice, as those considered in this work. 
Therefore, our first task, is to generalize this technique to include this class of systems. In order to 
accomplish this goal, we must face two problems: on one hand, due to the non-Bravais character 
of the lattices, we cannot characterize the full system with only one order parameter; on the other 
hand, the finite clusters considered have to have a minimum size such that all the order parameters 
variables are included in the reduced trace, at least for one of the clusters, which involves a more 
complicated mathematical problem. These difficulties, of course, are not a particular flaw of the 
EFRG scheme, as they are present in all the real space RG schemes, and are related to the fact that 
non-Bravais lattices are non-invariant under a scale transformation. To cite an example, under a 
block transformation, the kagome lattice does not transform in itself, but in a triangular lattice. 

The easiest way of solving these technical difficulties is by subdividing the non-Bravais lattice 
in a number of Bravais lattices, and allowing for a different order parameter in each sublattice. For 
example, in the case of the kagome lattice, it can be considered as formed by n = 3 interlocked 
triangular lattices (a = A, B, C), whereas in the pyrochlore case, it can considered as formed by 
n = 4 interlocked FCC lattices (« = A, B, C, D). 

Next, we have to decide the size of the two finite clusters we are considering. The most obvious 
choice is to consider a cluster formed by one single spin and another one formed by p = 3 spins 
for the kagome lattice (see Fig. and p = A spins for the pyrochlore lattice (see Fig. 

The following step is to evaluate the expectation value of each of the order parameters for these 
clusters by using the Callen-Suzuki identity, a task that we carry out in the next section. 



Let us consider the cluster defined by the ion 1 in sublattice a. The Hamiltonian of this cluster 
is given by 



III. EVALUATION OF THE ORDER PARAMETERS 



A. 1-spin cluster order parameter 



z 




(16) 



/37^a i 



8 



where the subindex la makes reference to the fact that the Hamikonian corresponds to a cluster 
of 1 spin belonging to sublattice a, and stands for the symmetry breaking field acting on the 
spin of the 1-spin cluster, which belong to sublattice a. The sum over the i index is performed 
over the NN of the la spin, with z the number of NN in each of the sublattices (z = 2 for both the 
kagome and pyrochlore lattices). By using the Callen-Suzuki identity we can calculate the order 
parameter in sublattice a for a cluster formed by 1 spin 



mia = {Sia) 



(17) 



H 



As commented above, Tria represents the partial trace with respect to the variables of the spin 
la, and the symbol (• • • )^ stands for the usual canonical thermal average taken over the ensemble 
defined by the complete Hamiltonian H. This equation can be put in the more convenient form 



m 



la 



d 



■ In TriaC 



-Hi 



where stands for the /x-th component of the vector a. A simple calculation yields 



(18) 



Next, we make use of the operator identityEl 

f{x + a) = 



fix) 



where = -§^, to rewrite eq. (|T9|) in the more convenient form 



mio, = Tic 



^n^-iio,D^\ ^anhxl 



n 



x=o ■ 



Finally, by making use of the van der Waerden identity 

e"'^ = cosh a + S sinh a, S = ±1, 



(19) 



(20) 



(21) 



(22) 



we can put 



exp 



Sipl'x 



(23) 



where = cosh (JC'Dj.) and = sinh {IC'D^), with fC' = K' cos9 = K'n^ ■ n^. The expression 
of K,' takes into account the fact that, for the structures considered in this work, the angle between 
the directions of anisotropy is the same for any pair of neighboring spins. 
Then, the order parameter for the 1-spin cluster can be put in the form 



mia = na{ n n [Pi 



tanhxl 



(24) 



H 
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For example, for the kagome lattice, the previous expression takes the form 

{a ^(3^^). (25) 

These kind of equations, though exact, are intractable, and we have to resort to some kind of 
approximation. The simplest one corresponds to neglect 2-spin and higher order correlations^^, 
that is 

{SiaSjfS ■ ■ ■ Sfc^)^ — (SjQ,)^ ('^i/S)^^ • • • ('^'27)7-^ ) (26) 

which is equivalent to neglect terms of the order of K"^ and higher in the calculation of the critical 
points. It must be noted, however, that this approximation is superior to the standard MF theory, 
since it takes exactly into account the relation (sf^) = 1 (i ^ 1), through the van der Waerden 
identity. 

Applying to (P^, and to first order in (sia)-^ = Wa (which constitutes a good approxi- 
mation near the critical point), where the subindex 1 means that these symmetry breaking fields 
have been obtained from the 1-spin cluster, we arrive to the final form of the order parameter in 
sublattice a 

mice = ria A{IC')J2bip, (27) 

where 

A{}C') = zp^(f-i)-V, tanhxl^^o , (28) 
which, after a straightforward calculation it is found to be 

^ 2 tanh(2/C0 + tanh(4A:0 ^ ^^^^ 

for the kagome lattice (z = 2 and p = 3), and 

^ 5 tanh(2/C') + 4 tanh(4/C') + tanh(6/C') 

^(^ ) = ' (30) 

lb 

for the pyrochlore lattice (z = 2 and p = 4). 



B. p-spin cluster order parameter 

Let us consider now a cluster formed by p spins, each of one belonging to a different sublattice 
(therefore, we will have p = 3 for the kagome lattice and p = 4 for the pyrochlore lattice). The 
Hamiltonian of such a cluster can be put as 

Hp = K^^Sia- Sip + ^ Sia ■ ipa- (31) 
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In this Hamiltonian, the first term represents the interaction between the spins inside the cluster of 
size p, whereas the second term represents the interaction of those same spins with the symmetry 
breaking fields 

2-1 

^p^ = kY,^s,p. (32) 

P^a i 

created by the rest of spins outside the cluster. The sum over i goes over all the NN spins except 
the ones already included in the cluster. The order parameter for each of the sublattices included 
in the cluster, obtained from the p-spin cluster can be now put as 

/TriA,iB,ic...Siae'^''\ 
mpa = i -7^ ^ ) , (33) 

or, more conveniently 

^^:a = (7JrlnTri^,iB,ic...e^A . (34) 
\ '^spa / n 

In these expressions, Tr\A,\B,\c... represents the trace over the spin degrees of freedom in the 
p-spin cluster. 

Even though the discussion can be continued for a cluster with an arbitrary number of spins, the 
mathematical details get very cumbersome. Therefore, we will carry on separately the calculation 
of the order parameter for the kagome and the pyrochlore lattices. 



1. Kagome lattice 

In the case of the kagome lattice, we will consider a triangular cluster defined by p = 3 spins. 
In doing so, we can easily calculate the reduced partition function appearing in (R^ 



TnAABAce'^" = ZsiJC, Ua, Ub, Uc) = cosh(n^ + Ub + R 



c) 



+ e"^^ [cosh(nA - + He) + cosh(n^ + - Uc) + cosh(nA - - He)] , (35) 
where IIq = ■ Csa and fC = K cosh 6. Taking into account the relation 

= <, (36) 

we obtain the expression of the order parameter for the 3-spin cluster 

' sinh(n^ + Hb + Tic) + e-^^ [± sinh(nA - + Tic) 



ZsilC,UA,TiB,Tic) 

± sinh(n^ + Ub- Tic) ± sinh(n^ - - n^)] \ ^ 
+ Z,{IC,UA,UB,Tic) In""' ^ ^ 
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where the + or — sign corresponds to the sign of IIq, in the argument of the corresponding sinh. 
To quote an example, the order parameter in sublattice B is given by 



'm3B 



sinh(nA + Ub + Uc) + e-^'^ [- sinh(nA -Ub + R 



c 



+ 



Z3{}C,Ua,Ub,Uc) 
+ sinh(nA + IIb- Uc) - sinh(nA - - Uc)] 



UB, (38) 



Zs{}C, Ua, Hb, He) / ^ 

By making use again of the operator identity (|20|), we can express the order parameters as 



.UaD^ UBDy UcD, 



)n 9c{x,y,z)\ 



x=y=z=0 ' 



(39) 



where ga{x, y, z) is obtained by making the substitutions 11^ — x, — > and Tic z in (|T7|). 
For example 



gB{x,y,z) = 

sinh(x + y + z) + e~^^ [— sinh(x — y + z) + sinh(x + y — z) ~ sinh(x — y — z)] 
cosh(x + y + z) + e""''^ [cosh(x — y + z) + cosh(a; + y — z) + cosh(x — y — z)] 



(40) 



By expanding the exponentials of sums in ( |39D as products of exponentials, and making use 
again of the van der Waerden identity and the decoupling of the correlations, we arrive, to first 
order in the symmetry breaking fields b^^ = (■SiQ)^^^ to the following system of equations 



'37 ) J IT'a ) 



("7^/5 7^7), 



where 



and 



^l(^) = plpyPziVyPz + PyVz) gA{x,y,z)\ 



x=y=z=0 ' 



52(/C) = plpxpzii'xpz + Pxi'z) gA{x,y, z) 



\x=y=z=0 ' 



(41) 



(42) 



(43) 



with = cosh{ICD^) and z/^ = sinh(/CD^), = x,y, z. After applying the differential operators 
in the previous expressions, we arrive to the final form for these coefficients 



5i(/C) 



1 
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1 



2 cosh(4/C) + Se-^'C 



+ 



6e 



cosh(4/C) + e-4'=(2 + cosh(4/C)) 

-4/C 



l + e-4'=(l + 2cosh(4/C)) 



sinh (4/C) + 



4tanh(2/C) 



(44) 



+ 



3(2 + e- 



2 cosh(4/C) + 3e-4^ - 1 cosh(4/C) + 6-4^^(2 + cosh(4/C)) 



3e 



1 + 6-4^(1 + 2 cosh(4/C)) 



sinh(4/C) + 



4(2 + 3e~^'^) tanh(2/C) 
1 + 3e-4^ ' 



(45) 



Certainly, it is not obvious to see the behavior of Bi and B2 from the above expressions, for they 
are depicted in Fig. 0. 
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FIG. 4: /C dependence of BiQC) (solid line), B2QC) (long dashed line), Ci(/C) (short dashed line), and 



C2(/C) (dotted line), given by (H3), (0|), (49), and dsq), respectively. 



2. Pyrochlore lattice. 

The only difference with the previous case is that the cluster contains p = A spins, correspond- 
ing to 4 different sublattices {A,B,C, and D). Therefore, we will not present the details of the 
calculation here, but only the main results. 

The reduced partition function is given in this case by 

Z4(/C, Ha, Hb, He, Hd) = cosh(nA + IIb + ^c + ^d) 

+ e"^^ [cosh(nA - + + Hz)) + cosh(nA + YiB-Tic + Hz?) 



+ cosh(nA + + He - Hfl) + cosh(nA - Hb - He - n^)] 

rcosh(nA - - + Hd) + cosh(nA - + - n^,) 

+ cosh(nA + Hb - He - n,,)] 



where the definition of IIq is as given above. 

The expression of the order parameter for this cluster is given by 



sinh(nA + Hb + He + Yip) 

Z4(/c,nA,nB,nc,nB) 

/± sinh(nA - Hb + He + n^^) ± sinh(nA + Hb - He + Ho) 



+ e- 



Z4(/C, n^, Hb, He, Hz,) 
± sinh(nA + Hb + He - Hb) ± sinh(nA - Hb - - Hg) 

Z4(/c,nA,nB,nc,nB) 

.gyc /± sinh(nA - Hb - He + Hb) ± sinh(nA - Hb + He - Hb) 



Z4(/c,nA,nB,nc,nB) 

±sinh(nA + nB-nc 



n 



Z4(/c,nA,nB,nc,nB) 



(46) 



(47) 
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where, as before, the + or — sign corresponds to the sign of the corresponding IIq, in the argument 
of the sinh. 

By making use of the operator identity <^^, the van der Waerden identity, and the decoupling 
approximation, we arrive to the expression of the 4-spin cluster order parameter up to first order 
in the symmetry breaking fields 64Q, = (sjQ,)^, which is given by the system of equations 

m4a = [C^i(/C) 64a + C^2(/C) (64/3 + ^47 + &45)]n„, {a^(3^-f^6), (48) 

where 

Cl{IC) = plplplpl{Vypzpu + PyVzpu + Pypz^u) QaIx, V, Z, u)\^^y^^^^^Q , (49) 

and 

C2(^) = pIpIpIpH^^xPzPu + Px^zpu + PxPzVu) 9a{x, y, z, u)\^^y^^^^^Q . (50) 

The function gaix, y, z, u) is obtained from the coefficient of Ua in ( p7[ ) with the substitutions 
Ha — ^ X, YIb — ^ y. He — > z, and YId u. 

The evaluation of Ci and C2 is straightforward, though rather lengthy, and the final expressions 
are of little use, due to their length. For that reason, we have preferred not to quote them here, but 
illustrate their behavior by Fig. ^. 



IV. SCALING RELATIONS AND CRITICAL BEHAVIOR 



Let us summarize the results obtained so far. The order parameters of each sublattice for the 
1-spin cluster are related by the system of equations 



(51) 



whereas the order parameters of each sublattice for the p-spin cluster are related by the equations 



m. 



pa 



$(/C)6p« + e(/C)5^6p^ 



rir 



(52) 



where A(/C'), $(/C), and e(/C), are given by expressions (g9|), (Q, (g5|), and (g|), (|g), for 
the kagome and pyrochlore lattices, respectively. 

By making use of the scaling hypothesis (Q) for both the order parameter and the symmetry 
breaking fields§ we arrive to a system of equations for the symmetry breaking fields 



$(;C) bp^ + [e(/C) - A{JC')] J2 bpP = 0, (53) 

where $, 6, and A are to be replaced by the corresponding expressions, and a, P = A, B, and C 
for the kagome lattice and a, P = A, B, C, and D for the pyrochlore lattice. This homogeneous 
system has a non trivial solution if its characteristic determinant is zero, which gives the condition 



[{p - 1)A(/C') - $(/C) - - 1)0(/C)] [A{}C') + <I>(/C) - Q{IC)r' = 0, (54) 
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(a) Kagome lattice. 



(b) Pyrochlore lattice. 



FIG. 5: Type A configuration for the kagome and pyrochlore lattices at T = 0. The ^ = case corresponds 
to the standard Ising model, whereas 9^0 corresponds to anisotropic Ising interactions (see text). 



which has two solutions that give us the corresponding RG recursion relations for the coupling 
constants. 

The first one corresponds to 



<I>(/C) + ip- l)e(/C) = (p - l)A{K') 



(55) 



In this case, the solution of the system of equations is given by (notice that we have omitted the 
index corresponding to the size of the cluster) 



bA = bB = bc = ... 
or, in terms of the modulus of the order parameter, 

rriA = mB = mc = . . . , 



(56) 



(57) 



and corresponds to a configuration in which all the spins point "out" or all the spins point "in", as 
has been depicted in Fig. |5|. We will term this configuration "Type A". 
The second solution occurs if 



G(X:) - <I>(/C) = A{}C') 
which leads to a configuration in which 



(58) 



(59) 



or, in terms of the order parameters 



= 0. 



(60) 



We will term this configuration as "Type B". 

It is important to notice that through the introduction of the effective couplings JC' and /C, 
we can consider both the standard Ising Hamiltonian in these geometries and the corresponding 
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(a) Kagome lattice. (b) Pyrochlore lattice. 

FIG. 6: Type B configuration for tlie kagome and pyrochlore lattices at T = 0. The = case corresponds 
to the standard antiferromagnetic Ising model, whereas 0/0 corresponds to anisotropic Ising interactions 
(see text). 



spin ice problem on equal footing or, equivalently, we can map the spin ice problem in these 
geometries to the standard Ising model. Indeed, the standard Ising model corresponds to taking 
cos 9 = na-np = I, that is, the spins point along the same direction, whereas the case 9 = 120° for 
the kagome lattice corresponds to the 2D spin ice analog commented above. The case 9 = 109.47° 
is the spin ice problem in the pyrochlore lattice. 

In the case = 0°, the Type A configuration corresponds to ferromagnetic order, whereas the 
Type B configuration is a generalization of antiferromagnetic order. Moreover, in the pyrochlore 
case, and for Ising spins pointing along (1, 1, l)-type directions, the last configuration contains, 
among others more general, the ground state configuration of the spin ice system, in which two 
spins point "in" and two spins point "out" (see Fig. H), and reflects the degeneracy of the ground 
state due to the geometrical frustration. It is important to notice that it impossible to satisfy condi- 
tion ( |60| ) in the kagome case at T = 0, which has been noted in Fig. |§ by a double spin in one of 
the comers. 

Of course, apart from obtaining the possible spin configurations, it is even more important to 
study if they are thermally stable, that is, if they really occur. To this end, we have to see if there 



is any value /C' = /C = /Cc that satisfies conditions ( p5| ) and (p8|), that is, if the RG recursion 
relations defined by ( |55| ) and ( |5^ ) have any non-trivial fixed points. The simplest way of solving 
these equations is by using graphical methods, as has been done in Figs. ^ and |S| for the kagome 
and pyrochlore lattices, respectively. 

From these figures, we can see that for the Type B configuration, the only fixed point of the RG 
transformation is the trivial one, corresponding to a paramagnetic configuration (/C = /C' = 0). 
However, there is a non trivial fixed point for the Type A configuration (/Cc = 0.415 for the kagome 
lattice and K-c = 0.23 for the pyrochlore lattice), depending on the sign of the effective coupling 
or, in other words, depending on whether the effective coupling is ferro- or antiferromagnetic. 
These results have been collected in Table |, for both the kagome and pyrochlore lattices, and for 
both the standard Ising model and the spin ice problem. Also, the exact results known for the 
kagome latticelii and the ones obtained from MC data for the pyrochlore lattice! are quoted for 
comparison. It can be seen that the estimates obtained from this RG scheme are in very good 
agreement with the exact values, in spite of the crudeness of the decoupling approximation and 
the fact that we have considered clusters of small size. This reflects the well known fact that, in 
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(a) Type A configuration. 
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(b) Type B configuration. 



FIG. 7: Graphical estimation of the fixed points of the RG transformation for the kagome lattice. A{x) is 



given by (29), whereas Bi{x) and B2{x) are given by (42) and (43), respectively. 



JCc Kagome Pyrochlore 
EFRG 0.415 0.230 
GCC 0.402 0.227 
Exact 0.4643'' 0.25^ 



''Ref. 
''Ref. 



TABLE 1: Non-trivial fixed points of the RG transformation, and corresponding critical points in the GCC 
approach (see text). 



this GFAF systems, correlations extending outside the frustrated unit are relatively sma ll§lil. On 
the other hand, correlations inside the unit are considered to some extent by the present technique. 
In any case, the quality of the estimates can be systematically improved by considering larger 
units as, in this way, we include further correlations. However, the improvement in the estimate 
of the critical point is not commensurate with the additional algebraic difficulty involved. Another 
important point is that the EFRG method distinguishes among the various different geometries for 
a given number of NN, in contrast with standard MFapproaches, which do not distinguish, for 
example between a cubic lattice and a pyrochlore oneH3. 

Regarding the question of geometrical frustration, which leads to the absence of a transition to 
a long range ordered state at any finite temperature, both the kagome and pyrochlore ferromagnetic 
Ising lattices are non-frustrated, whereas the corresponding antiferromagnetic ones are frustrated, 
and remain paramagnetic down to K. In the spin ice case (9 = 120° for the kagome lattice 
and 9 = 109.47° for the pyrochlore), however, the ferromagnetic case is frustrated, and does not 
order at any finite temperature, whereas the antiferromagnetic one experiments a transition to a 
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Kagome 


Pyrochlore 






A > U A < U 


A > -K < 


EFRG 


^ = 


0.415 — 


0.230 — 




e = 120°(109.47°) 


— -0.830 


— -0.690 


GCC 




0.402 — 


0.227 — 




e = 120°(109.47°) 


— -0.804 


— -0.671 



TABLE II: Critical points obtained from the EFRG and GCC methods for the Ising, kagome spin ice, and 
spin ice pyrochlore. The long dashes indicate the absence of a critical point. 



1 

0.8 
0.6 
0.4 
0.2 


-0.2 
-0.4 
-0.6 
-0.8 



































-Gi(x)/3+ 








G2(x) 































-0.5 





X 



0.5 



0.8 
0.6 
0.4 
0.2 


-0.2 
-0.4 
-0.6 
-0.8 



















































































C2(> 
































A(x 


) 















(a) Type A configuration. 
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(b) Type B configuration. 



FIG. 8: Graphical estimation of the fixed points of the RG transformation for the pyrochlore lattice. A{x) 
is given by (30), whereas C\{x) and C2(x) are given by (p^) and (50), respectively. 



long range ordered state at a finite temperature (see Table In any case, the ground state given 
by the "ice rule"LKl never occurs in these systems. Furthermore, we can see now that the reason 
why the ferromagnetic spin ice system is frustrated is that it can be mapped to the standard Ising 
Hamiltonian with antiferromagnetic effective interaction, /C < 0, which is known to be frustrated. 



V. CALCULATION OF THE CRITICAL EXPONENTS 

In this section we will carry out the calculation of the critical exponents associated with the 
transition to the Type A configuration, which, as found in the previous section, is the only one that 
occurs in these systems. As it is easy to see from the considerations in Section ||, in order to obtain 
these exponents, we will need to include the effect of a magnetic field applied to the system, which 
is what we do in the following. 
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We can make use of the results of the previous sections to consider the inclusion of the applied 
magnetic field, h = hn, in the direction given by the normal vector n, by redefining the conjugated 
fields as 

— * — # 

^pa ^pa + {n ■ Ua) h (61) 

for the p-spin cluster and 

iia fiQ + {n-na)h' (62) 
for the 1-spin cluster. The magnetic fields for both clusters scale as 

h' = ly^h, (63) 

with uh the magnetic critical exponent0. 

Now, it is straightforward to see that the corresponding order parameter in each sublattice for 
the 1-spin cluster is given by 



fhia = ria ^tanh 



na ■ ( 6a + nh 



n 



For small values of h, we can put 



mia — "n-Q, <5 ^tanh 



+ ^sech^ 



la 



{n-na)h + 0{h'^ 



(64) 



(65) 



By making use of the exponential operator identity together with the van der Waerden identity, 
we easily obtain, to the lowest order in the fields and h 



mia 



(66) 



where A{fC') is given by ( P8| ) and 



F(/C') =p'}'P-^^ sech^l 



X) 



x=0 



(67) 



For the kagome lattice (z = 2 and p = 3) 



F(/C') = -(3 + 4sech2 2/C' + sech2 4/C'). 



(68) 



For the pyrochlore lattice (z = 2 and p = 4) 



F(/C) = — (10 + 15sech^2X: + 6sech2 4/C + sech26/C). 
32 



(69) 



Let us now turn our attention to the 3-spin cluster in the kagome case. We follow exactly the same 
procedure we have used to arrive to eq. (|5^, but using the expression for the order parameter of 
this cluster, eq. (pT]). Thus, we first expand rh^a in powers of h and b^a, after the substitution ( |^ ) 
has been done, and make use of the differential operator and van der Waerden identities. In doing 
so, we arrive at 



msaihA, bsB, he, h) ~ rh^a{hA, hs, he) + na{Ga{10j ■ h) + 0{h'^ 



(70) 
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The first term on the right member of this expression is the same as (|4T|), whereas the second one 
is new, and the function Gq,(/C) is given by 

d 



dh 



{x + ua ■ h, y + ub ■ h, z + fie ■ h) 



(71) 



h=0 



with ga{x, y, z) given as in and (PPj). For arbitrary n^, tib, and ric, Ga will depend on a in 
a complicated way. However, in this work, we are only concerned about the cases studied in the 
previous Section, namely, the standard Ising model, ua = n-B = nc, and the 2D spin ice problem 
on the kagome lattice, which verifies the property nA + ns + nc = 0. In these cases, it is easy to 
show that 



where the form of the function G{}C) does depend on whether ua 
In the first case (standard Ising model), this function is given by 



GiilC) 



192 

-6 + e^^) 



138- 



320 



+ 



704 



6144 



;i + e4'C)2 i + e^'c (3 + e4'C)2 
3(17 +lle^'^) 3(1 + 806^^^) 



+ 



Ub 



768 



(72) 

nc or HA+nB+nc = 0. 



3 + e 



4/C 



+ 



24(-7 + 3e^'^) 
(4 - e^'^ + e 



8/C^2 



4 - e^'c + (1 + e^^ + 2e8'C)2 



_ 48(-13 - 50e^^ + 636^'^) 
1 + e^^ + 2e8^ ~ (1 + Se^'c + + e^'^^f 
48(-13 + 18e^'^ + e^^; 



+ 



1 + Se^'^ + eS'c + e 



\2K 



(73) 



whereas, in the second case (2D spin ice), it is given by 

GsiijC) = sech^ 2/C (533575 + 1162616 cosh 4/C + 1014592 cosh8/C + 774088 cosh 12/C 
+436732 cosh 16/C + 178200 cosh 20/C + 47360 cosh 24/C + 6824 cosh 28/C + 317 cosh 32/C) / 
[1024(3 cosh 2/C + cosh 6/C - 2 sinh 21Cf{-2 cosh 2/C + sinh 21Cf 

X (-1 + 5 cosh 4/C - 3 sinh 4/C) (1 + 3 cosh 4/C + sinh AlCf] 
- sech^ 2/C (440942 sinh 4/C + 661942 sinh8/C + 618862 sinh 12/C 
+385138 sinh 16/C + 163174 sinh 20/C + 44430 sinh 24/C + 6374 sinh 28/C + 303 sinh 32/C) / 
[1024(3 cosh 2/C + cosh 6/C - 2 sinh 2/C) ^(-2 cosh 2/C + sinh 21Cf 

x(-l + 5cosh4/C -3sinh4/C)(l + 3cosh4/C + sinh4/C)2] . (74) 



The /C dependence of these expressions has been depicted in Fig. ^ 

The calculation for the pyrochlore lattice is identical to the one of the kagome, though the 
algebraic details are a little more involved. The final result of the calculation is qualitatively the 
same as in the kagome case considered above: the magnetization of the 4-spin cluster can be put 
as 

m4a(64A, hiB, he, ho, h) ~ m4„(64A, &4B, he, ho) + H{1C) (fia ■h)fia + 0{h^). (75) 



where the first term on the right member is given by ( PS| ) and H{1C) is obtained from the ga 
function for the pyrochlore case as 

HailOj = pIpIpI —ga{x + nA- h,y + nB ■ h,z + nc ■ h) 



(76) 

h=0 
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FIG. 9: K, dependence of the functions Gj{IC) (solid line), Gsi{IC) (long dashed line) (given by ( [73| ) and 
(0), respectively), Hj{IC) (short dashed line), and Hsi{K,) (dotted line), (obtained from ([7^, see text). 



Again, for the cases considered in this work, ua = ub = nc = n-D (standard Ising) or ua + 
ub + nc + Uf) = (spin ice), it can be easily shown that Ha{)C) = H{K,) Ua, where the form of 
the function H{IC) does depend on whether we are considering the standard Ising or the spin ice 
problem. We will represent by Hi{K,) the first possibility and by Hsi{JC) the second. The length 
of these expressions makes them of little practical use and, therefore, we will not quote them here. 
It is enough for our purposes to give their /C dependence by means of Fig. |^. 

Now, by using the scaling relation (Q), we can equate the results for the 1-spin cluster with 
those for the p-spin cluster, as we did in Sec. In doing so, we obtained the following system of 
equations for the internal fields 

$(/C) + [e(/C) - A{IC')] J2 bpP = r^+'y^ - S(/C)] {K ■ n)h, (77) 



where $(/C), e(/C), and A{IC) are the same as in Sec. f^, F(/C') is given by and for the 
kagome and pyrochlore lattices, respectively, and S(/C) is given by G{1C) {Gj{IC) for the standard 
Ising model, and Gsi{K,) for the 2D spin ice) for the kagome lattice, and H{K) (if/(/C) for the 
standard Ising model, and Hsi{JC) for the spin ice) for the pyrochlore lattice. 

In the limit case /i = 0, we recover the results obtained in Section for the critical point. How- 
ever, now we have an additional relation that allows us to evaluate the magnetic critical exponent 
as 



(78) 



Moreover, the correlation length critical exponent is calculated by making use of ([T^) and ([13]) as 

In I 



In Xn 



(79) 
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Kagome Pyrochlore 





V 




V 


yu 


Standard Ising (ua = nB = 


...) 1.07 


1.64 


0.83 


2.32 


Spin Ice {J2a = 0) 


1.07 


0.63 


0.83 


1.22 


Exact (Standard Ising)" 


1 


1.875 


0.63 


2.48 



Refs.||,|3ig|5| 

TABLE III: Values of the critical exponents. 



where 

~ o "I r "I — ^ 

\t= ^($(/C) + (p-1)0(/C)) (p-1)^A(/C') . (80) 



^ ($(/C) + (p-l)0(/C)) 



Sir 



The values of z/ and yu obtained for both standard Ising spins and spin ice spins, and for both 
the kagome and the pyrochlores lattices, are quoted in Table The other exponents for the 
different thermodynamic quantities can be obtained from these by making use of the well known 
scaling relationstj (see Table |[V|). The values listed as "exact" are those obtained for the square 
and cubic lattices, as they are expected to be the same as for the kagome and pyrochlore lattices, 
respectively, on the basis of the universality hypothesis. The ones corresponding to the square 
lattice are those obtained from the Onsager solutionE3 wdiereas the ones for the cubic lattice are 
values obtained from high-precision numerical studiesc3tSB. 

The first conclusion we can extract from the listed values for the standard Ising problem is 
that the EFRG method provides very reliable estimates for the critical exponents, in spite of the 
fact that we have chosen the smallest possible clusters that include the effects of frustration. In 
this sense, the exponents calculated in this work are some of the "best" that can be found in the 
literature, calculated in the framework of phenomenological RG approaches. For example, the 
exponents obtained in our work are even better than the ones obtained in the so called Mean Field 
RG method for Ising spins in the square lattice, where clusters of 16 and 25 spins where usedEl. 
The reason for this improvement is in the use of the van der Waerden identity, which exactly 
takes into account the value of the autocorrelation function. It is also important to note that the 
estimation of v for the kagome lattice is better than the one for the pyrochlore lattice. The reason 
is that the scaling factor we have used, I = (N/N'Y^'^ = p^^'^, is smaller for the 3D lattice than for 
the 2D one. It has been noted by some authors, that the values of the critical exponents calculated 
with this method can be improved by using an alternative definition of the size of the clusters, in 
terms of the number of interactions including the ones with the spins creating the effective fieldE3. 
This point will be further investigated in a future workEJ. 

The last important conclusion we can extract from Table ^ is that the magnetic critical ex- 
ponent and, thus, all the exponents depending on this one (P, 7, S, and r]), are dijf event for the 
spin ice case. If we take into account the good agreement found for the standard Ising case, we 
can expect the estimations of the critical exponents for the spin ice case to be of the same quality. 
Actually, if we assume that the error in the determination of the critical exponents for the spin ice 
case is the same as in the standard Ising case {/^yulVH = 0.13 for the kagome Ising case and 
^VhIvh = 0.06 for the pyrochlore Ising), and we see no reason for this not to be the case, we 
can put an upper and lower bounds for the real value of this exponent, yu = 0.63 ± 0.08 for the 
kagome spin ice, and yn = 1-22 ± 0.07 for the pyrochlore spin ice. Apart from a possible relative 
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^ - d-IVu 


R 
P 


{d — yH)/yT 


7 


{2yH - d)/yT 


5 


VhI {d - yu) 


V 


2-2yH + d 


V 


l/Vt 



TABLE IV: Expressions of the critical exponents in terms of the RG critical exponents and yn- 



error of about a 10% in the worse case, they are obviously very different from the of the stan- 
dard Ising model and, in fact, we have not found in the literature any reference to any other system 
with these values for yu- A naive interpretation of these exponents could make us think that our 
results violate the universality hypothesis. However, it is easy to realize that this is not the case. 
Indeed, the universality hypothesis states that the critical exponents are the same for different sys- 
tems with Hamiltonians which internal symmetries are the same, independently of details as the 
microscopic origin of the interactions or the geometry of the lattice. The results we found for the 
standard Ising model in both the kagome and pyrochlore lattices are in good agreement with this 
idea. Moreover, the fact that the thermal exponent does not change when going from Ising spins 
with a single anisotropy axis, to the spin ice problem, is also in agreement with the universality 
hypothesis, as we have previously shown that the spin ice problem maps on an Ising Hamiltonian 
with effective interactions given hy }C = K cos 9. However, when we include the effect of an 
uniform magnetic field in the spin ice case, the Zeeman term cannot be mapped to the Zeeman 
term in the Ising Hamiltonian with a single anisotropy axis, which is the same for all the spins, 
but we can see it as a Zeeman which is different for spins located in different sublattices, i.e., it 
is equivalent to some kind of staggered field subject to the condition J2a = 0, where by ha 
we represent the magnetic field acting on each sublattice. Certainly, this Zeeman term has not the 
same symmetry as the original Zeeman term of the Ising spin model with a single anisotropy axis 
and, therefore, we should expect the corresponding magnetic critical exponent to be different from 
the one of the standard Ising model. Of course, it would be extremely important to have exper- 
imental confirmation of this fact, or even estimations of the magnetic critical exponent obtained 
from high precision Monte Carlo simulations but, as far as we know, this is the first time that this 
result is reported. 

Another interesting point we waiitip notice here is the connection of these spin ice-like systems 
with the well known vertex modelsBH. In fact, it is easy to see that the 2D spin ice of this work 
is related to the 8-vertex model, whereas the spin ice on the pyrochlore is related to the 16-vertex 
model. On the other hand, it is well known that these vertex models have the peculiarity that the 
critical exponents depend continuously on the value of the exchange coupling, and this is precisely 
what we would expect in our model if we allow for arbitrary orientations of the unit vectors n^, 
for which ua ■ riB ua -ub 7^ .... We think that clearly establishing this connection, though 
outside of the scope of this work, deserves further consideration, as it seems that these spin ice-like 
systems are one of the few physical realizations of these vertex models0. In this sense, we think 
that the present RG scheme could prove to be useful in obtaining qualitative information about 
critical properties of the 16-vertex model, about which very little is known. 
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VI. THE GCC METHOD APPLIED TO THE SPIN ICE PROBLEM 

As stated in the Introduction, it is interesting to compare the results obtained in the previous 
sections with those obtained in the framework of the GCC method, which has been successfully 
applied to the problem of Heisenberg frustrated lattices by the present authorsBEl. In the GCC 
framework, we initially consider isolated units (tetrahedra or triangles, for the pyrochlore and 
kagome lattices, respectively), and later add the interactions with the surrounding units as unknown 
effective fields to be determined by a self consistent conditioiL^hich is obtained by comparing 
the order parameter in the cluster with that of the isolated ionoEZl. We will analyze here the case 
in which the applied magnetic field is zero, even though the generalization is straightforward. 

As in the previous sections, we will characterize each spin in the frustrated unit by a different 
order parameter. Therefore, the Hamiltonian of the cluster formed by one Ising spin is given by 

Hla = 2K Sia ^hp = 2Sia " ^la, (81) 

where are the internal fields created by the neighboring spins outside the unit, and point along 
the directions given by n^, and we have taken into account that there are 2 neighbors in each 
sublattice. The Hamiltonian of the cluster formed by p spins can be put as 

Hp = K^Sia ■ Sip + ^Sia- ipa- (82) 

a^l3 a 

We can then see that the corresponding partition function of the unit has the same functional form 
as eqs. (^7]) and (|47|) for the kagome and pyrochlore lattices, respectively, except for the fact that 
the average with respect to the full Hamiltonian is absent, and the conjugated fields appearing 
in Ua depend on the internal fields hp, rather than on the spin variables of the surrounding spins 
outside the cluster. As stated above, these internal fields are evaluated by imposing the self con- 
sistent condition of equating the magnetization per spin in the internal field with that of a unit in 
the internal field, which can be mathematically stated as 

(jlA,hB,hc,- ■ - ^ ■ (83) 

Obviously, this equation can only be solved numerically in the general case. However, near the 
critical point, we expect the internal fields to be very small, so we can expand the order parameters 
to first order in the internal fields. In doing so, we obtain 

mia^2lcY,hf3, (84) 



mp„ = 2/C [ (/C) /i„ + ^2 (/C) 5^ /i J , (85) 



where )C = K cos 6, as before, and 

l-e- 



A(/c) 



l+3e-*'C 

3(l_e-8K) 

l+4e-6K+3e-8K 
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and 



M^) = { 'sl.t-'^+e-B'c , (86) 



6k;_i_q„-8k: 



l+4e-b'^+3e 

where the upper (lower) expression corresponds to the kagome (pyrochlore) lattice. 
Then, we arrive to the following system of equations for the effective fields 

^IC)K + e{IC){h^ + h^ + ...) = 0, (87) 

where 

$(;C)=Ai(/C), (88) 

and 

' 1 



e(/c) = ^(/c) - ; = i4ra.-5e-B^ • (89) 

Again, the upper (lower) expression corresponds to the kagome (pyrochlore) lattice. 
This system of equations has two non-trivial solutions for 

<I>(/C) = 0(/C) (90) 

and 

«D(/C) = -(p-l)0(/C). (91) 

The first one corresponds to a configuration in which J2a = 0, whereas the second one corre- 
sponds to a configuration in which Ha = = he = ■ ■ ., that is, they are the same configurations 
we found in the EFRG formulation. It is easy to verify that the first kind of solution never occurs 
for any finite value of /C, whereas the second one occurs for a value of ICc = 0.402 for the kagome 
lattice and JCc = 0.227 for the pyrochlore lattice, in very good agreement with our previous results 
(see Table I). 



VII. CONNECTION BETWEEN THE GCC AND EFRG SCHEMES 

From the previous section, it seems clear that both approaches presented in this work are some- 
how related. In this section we will furnish this connection. We will consider a system with only 
one vectorial order parameter. The generalization to an arbitrary number of order parameters being 
straightforward. 

The Hamiltonian of the one spin cluster with only NN interactions can be put as 

n = si -fi = si • [^i^^ ' (92) 

where z is the number of NN. The corresponding order parameter can be obtained then from the 
Callen-Suzuki identity and, in general, it will be a function of the conjugated field 

mi = (^lnZi{^i)) =(%0)^=($^i^Ji^X^sl) , (93) 



n 
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where define the coordinate system. Obviously, we can expand the function F in powers of Si 
as 



mi 



dF^ 

dsi 



Si=0 



d'F„ 



2! dsi ■ dsj 



Si ' S ^ 



n 



En 



Si=0 



d'F,„ 



■ (9sj 



If we now apply the decoupling approximation to this expression 



(94) 



(95) 



where hi = (s?)^, and take into account that the variables in the derivatives are dummy variables, 
that is, we can put 



dF^is,) 



dsi 



Si=0 



the series can be resumed in terms of the effective internal field hi to give 

fhi ~ F{z K hi). 

Obviously, we can do the same with the p cluster order parameter, and the result is 

fhp ~ G{{z — 1) Khp), 



(96) 



(97) 



(98) 



where hp = {si)^ is the effective internal field obtained from the p-spin cluster. We can then 
easily see that the GCC self consistency condition corresponds to the particular case of the finite 
size scaling hypothesis in which one assumes 



and 



hpi 



hp. 



(99) 



(100) 



Thus, the GCC model corresponds to the lowest order of the EFRG method, in which correla- 
tions outside the cluster are neglected and the scaling relation is substituted by strict equality. 

Therefore, the value of the critical point obtained in both methods, to this order of approxima- 
tion, is the same, as it does not depend on the particular scaling coefficient. However, this has 
profound consequences for the calculation of other quantities, as the critical exponents, which do 
depend on the particular scaling relation, as have been noted by some authorsli^. The differences 
between the values listed in Table | are not related to the use of the EFRG or GCC methods, but 
to the fact that we have used the van der Waerden identity in the case of the EFRG. However, we 
could make use of this identity in the GCC approach, at least for Ising spins, and the results would 
be identical to this order of approximation (i.e. up to terms linear in Kc). 

It is also important to stress that the results of this section are not limited to the connection 
between the EFRG and GCC methods applied to GFAF, as we have not made any assumption 
about the sizes of the clusters. In this sense, this connection also holds for the original formula- 
tion of both the Constant Coupling and EFRG methods, in which clusters with 1 and 2 spins are 
compared^. In other words, the Constant Coupling method (in whatever formulation) is the low- 
est order approximation of the EFRG method, corresponding to neglecting multispin correlations 
and assuming that the anomalous dimension of the order parameter is zero. 
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VIII. CONCLUSIONS 

In this work we have studied the critical behavior of ferro- and antiferromagnetically coupled 
Ising spins with local anisotropy axes in geometrically frustrated geometries (kagome and py- 
rochlores lattices). The present problem, applied to the pyrochlore lattice with local anisotropy 
axes in high symmetry (1 1 l)-like directions corresponds to the spin ice problem. The critical 
points have been calculated by means of the effective field renormalization group technique, ade- 
quately generalized to deal with the non-Bravais character of the lattices, and by the generalized 
constant coupling approach. These methods allow us to evaluate both the critical points, if any, 
and the corresponding ground state configurations. The results obtained are in very good agree- 
ment with exact results known for the Ising Hamiltonian in the 2D kagome lattice and Monte Carlo 
calculations for the pyrochlore lattice. For ferromagnetic interactions, geometrical frustration pre- 
vents the formation of long range order, whereas for antiferromagnetic ones, there is a transition 
at a finite temperature to a long range ordered state in which all the spins of the unit point outside 
or inside the frustrated unit. These results can be traced back to the fact that this spin ice problem 
can be mapped to the standard Ising model with an effective interaction given hy JC = K cos 9, 
where 9 is the angle between any 2 anisotropy directions. Therefore, ferromagnetic interactions in 
the spin ice problem correspond to effective antiferromagnetic ones, and vice versa. Incidentally, 
the connection between the effective field renormalization group technique and the generalized 
constant coupling method is furnished. It is shown that the GCC method corresponds to the lowest 
approximation of the EFRG scheme, in which correlations outside the cluster are neglected and 
the anomalous dimension of the order parameter is taken to be zero. 

In addition, the critical exponents of these systems have been calculated for the transition to 
the long range ordered state that occurs for effective ferromagnetic interactions. The calculated 
values of the critical exponents for spins with parallel anisotropy axes (the standard Ising model), 
are found to be in very good agreement with both those obtained from exact results for the 2D 
square lattice and those calculated computationally for the cubic lattice, in spite of the small size 
of the clusters considered in this work. This fact is in agreement with the universality hypothesis. 
For the spin ice problem (in both the kagome and the pyrochlore lattice), the thermal exponent 
has the same value as for standard Ising spins, as both Hamiltonians can be identically mapped on 
to each other. However, the magnetic exponent is found to be different from that of the standard 
Ising model in the presence of an uniform magnetic field. This result is expected, as the Zeeman 
term in this case has different symmetry properties than in the standard Ising model. In any case, 
independent confirmation of this fact, by experimental or Monte Carlo simulations, would be 
desirable. 

We think that the techniques presented in this work will prove to be useful in studying other 
phenomena present in the spin ice system. For example, it would be interesting to study how the 
inclusion of long range dipolar interactions modify the present results, that is, if these additional 
interactions could remove the degeneracy of the ground state so a long range ordered state is 
formed at any finite temperature for the ferromagnetic case. So far, there are numerical results 
and mean field calculations which indicate that the degeneracy of the ground state is weakly lifted 
by long range dipolar interactions, and there is a finite temperature much lower than the energy 
scale defined bv the dipolar interactions at which a long range ordered state given by the ice 
rules is formedElEI, but we think it would be important also to support this result with analytical 
calculations that go beyond mean field theory. Also, the effect of a small amount of dilution in 
the crystalline lattice is an intriguing problem which can be dealt with in the present formalisms. 
Moreover, even though the extension is not trivial, the present EFRG scheme can be generalized to 
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deal with Heisenberg spins and, thus, shed some light in the fascinating phenomena of geometrical 
frustration. 

Of course, the EFRG scheme is not exempt of limitations. For example, the inclusion of fluc- 
tuating corrections coming from correlations further than NN interactions is not an easy task. 
However, there are some methods to deal with this problemE3, though they are highly technical 
and outside of the scope of the present work. This, and some of the problems mentioned in the 
above paragraph are the directions of our present research. 
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It is important to stress that this single axis Ising spin system cannot occur in a pyrochlore lattice, as it 
is incompatible with the cubic symmetry of the spatial group. However, it is a very instructive example, 
and will prove to be useful for testing the accuracy of the results presented in this work. 
-^2 This is the case even if we consider the smallest possible clusters. See, for example, I. P. Fittipaldi and D. 
F. de Alburquerque, J. Magn. Magn. Mat. 104-107, 236 (1992), where they studied the Ising model with 
bond dilution in both the simple square and kagome lattices in the EFRG framework by using clusters 
of 1 and 2 spins, respectively. 

Another possible definition of the scaling factor we can use consists in calculating from (79) the value 
of / that leads to the exact value of yt for the standard Ising model, and make use of this value to 
estimate the corresponding magnetic critical exponent. The values of / obtained in this way are / = 1.67 
and / = 1.42 for the kagome and pyrochlore lattices, respectively. The corresponding magnetic critical 
exponents are y// = 1.68 and yjj = 2.58, for the standard Ising model in the kagome and pyrochlore 
lattices, respectively, and yn = 0.60 and yjj = 1.13, for the 2D spin ice and the 3D spin ice, respectively. 
Even though this is only the case for very special cases, as the ones considered in this work, because 
spin ice systems with arbitrary unitary vectors are, in general, incompatible with the symmetries of the 
crystalline lattice. 



